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CALCULATIONS OF LAMINAR HEAT TRANSFER AROUND CYLINDERS OF ARBITRARY 
CROSS SECTION AND TRANSPIRATION-COOLED WALLS WITH 
APPLICATION TO TURBINE BLADE COOLING' 


By E. R. G. Eckert and J. N. B. Livincood 


SUMMARY 

An approximate 7nethod jor the development of flow and 
thermal boundary layers in the laminar region on cylinders 
with arbitrary cross section and transpiration-cooled walls is 
obtained by the use of Karman^s integrated momentum equation 
and an analogous heat-flow equation. Incompressible flow 
with constant property values throughout the boundary layer is 
assumed. The velocity and temperature profiles within the 
boundary layer are approximated by expressions composed of 
trigonometric functions. Shape parameters for these profiles 
and functions necessary for the solution of the bounaary-layer 
equations are presented as graphs so that the calculation for 
any specific case is reduced to the solution of two first-order 
differential equations. 

The method is applied to determine local heat-transfer 
coeflicients and surface temperatures in the laminar flow region 
of the transpiration-cooled turbine blades for a given coolant 
flow ratCj or to calculate the coolant flow distribution which is 
necessary in order to keep the blade temperature uniform along 
the surface. 

INTRODUCTION 

Transpiration cooling is a very effective means for keeping 
surfaces that are subject to a hot gas stream at a low tem- 
perature. For use of this method, the surface is fabricated 
from a porous material and a cooling fluid is blown through 
the pores. Along the outside surface the cooling fluid builds 
a film that insulates the wall from the hot gas stream. The 
transpiration-cooling method may ]>e applied to the cooling 
of structural parts in propulsion systems such as gas-turbine 
blades, combustion-chamber walls, and rocket nozzles. If 
a heated fluid is blown thi-ough the porous wall, the same 
method may be used to keep the surface temiierature of the 
wall at a value that is higher than the temperature in the 
outside flow. In this way, the method may be appliinl in 
de-icing of wings or other parts of airplanes and in the jiro- 
pulsion system. 

This report presents a method by which the heat tninsfer 
connected with transpiration cooling in two-dimensional 
laminar flow around bodies of arbitrary cross section can be 
calculated. It considers only the case where the fluid 
blown through the porous wall is the same as the one in the 
outside flow. Although the particular application con- 


sidered in this report is the cooling of turbine blades, the 
method itself may be applied to other applications as well. 
The procedure by which the calculation may be carried out 
for any particular application is described in an appendix 
of this report. 

The determination of the heat transfer is based on the 
calculation of the thermal boundary layer which builds up 
around any l)ody in a flow field. This thermal boundary 
layer is interconnected with the flow boundary layer for 
variable fluid properties dependent on temperature or super- 
imposed on the velocity boundary layer when the properties 
are independent of temperature. Tlie build-up of the flow 
boundary layer is determined by the pressure distribution 
around the body under consideration. For the type of 
])ressure distribution occurring in the aforementioned a])plica- 
tions, only approximate methods of calculation are direct 
enough for engineering purposes. For the determination of 
the thermal boundary layer in particular, two types of ap- 
proach are known. 

TJie first approacli was introduced by Kroujiline for the 
calculation of heat transfer on solid surfaces and presented 
in more detail by Frosling and others (ref. 1). In this meth- 
od, the flow boundary layer has to be knovvm before the ther- 
mal boundary layer can be calculated. 

Usually a method such as that introduced by von Karman 
(ref. 2) or Poldliausen (ref. 3), which fulfills the integrated 
momentum equation of the boundary layer, is applied for 
the calculation of the flow boundary layer. More recently, 
Wieghardt and Walz (ref. 4) have used, in addition to the 
momentum equation, an integrated energy equation, and 
Tetervin and Lin (ref. 5) have introduced a still more general 
integral condition for the boundary layer which may he used 
in such calculations. The use of these cx])ressions gives 
better agreement with exact calculations and with measure- 
ments in spec'ial cases, particularly in regions where the pres- 
sure increases in flow direction. Since, in tlie ap])lication 
considered, the regions of most concern are those where the 
pressure decreases, the integrated momentum equation, 
which is simplest to handle, will be used. Schlichting (ref. 
()) used this equation to calculate the flow boundary layer 
on a porous surface through which fluid is sucked in order 
to keep the boundary layer laminar or to prevent flow sep- 
aration. A paper l)y Dorodnitzyn (ref. 7) extended the 
method to include the effect of Mach number and of variable 
property values. However, this extension is developed only 
for zero heat transfer. The calculation of the flow boundary 


> Supersedes NACA RM E51F22, “Calculations of Laminar Heat Transfer Around Cylinders of Arbitrary Cross Section and Transpiration-Cooled Walls With Application to Turbine 
Blade Cooling,” by E. R. Q. Eckert and John N. B. Livingood, 1951. 
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layer in this re])ort will essentially follow Selilieliting’s 
approach. 

After the flow boundary layer is determined, the thermal 
boundary layer can be calculated according to the method of 
Kroujiline by use of an integrated heat-llow equation. Ve- 
locity profiles known from the calculation of the flow bound- 
aiy layer and temperature profiles within the boundary layer 
whose shapes are approximated by a polynomial expression 
are introduced into this heat-flow equation and the equation 
is solved for the thickness ratio of the thermal boundary 
layer to the flow boundary layer (ref. 1). This ])rocedure, 
however, becomes quite lengthy, especially when the number 
of terms in the polynomial expression for the tem])erature 
profile is increased in order to improve the accuracy of the 
ap])roxiniation. When an attempt was made at the N AC A 
Lewis laboratory to extend this method, which was previously 
used only to determine the heat transfer on solid surfaces, to 
transpiration-cooled ])orous walls, it was found that the 
procedure for the solution of the heat-flow equation becomes 
much simpler and more direct if the equation is solved for a 
thermal l)oundary-layer thickness termed convection thick- 
ness. This will be explained in detail in a later section of 
this report. After the method was developed and the cal- 
culations were finished, a note by Dienemann (ref. 8) was 
found in which the same kind of approach is briefly described. 
Dienemann applies the method to calculate heat-transfer 
coeflicients on solid surfaces and proposes to extend it in 
such a way as to account for a temperature variation along 
the solid surface and the influence of internal frictional heating 
within the boundary layer. He also shows that this method 
is superior to other apiiroximations which were compared 
l)y Goland (ref. 9) with an exact solution for a cylinder with 
circular cross section. It may be of interest to mention that 
Cioland obtained the exact solution from the fact that the 
dilferential equation describing the temperature boundary 
layer around an infinite cylinder in a flow normal to its axis 
and for a fluid with a Prandtl number of 1 has exactly the 
same form as the differential equation describing the span- 
wise flow within the boundary’ layer on a yawed infinite cyl- 
inder. Calculations of the heat transfer on a transpiration- 
cooled flat plate which included the variation of property 
values with temperature were made by Yuan (ref. 10) with 
the assumption that the total-temperature profile within the 
boundary layer is similar to the velocity profile. This as- 
sumption is valid for a Prandtl number of 1 and for no pressure 
gradient. The purpose of the present investigation is to 
consider the influence of ])ressure gradients as well as 
lYandtl numbers different from unity. 

A second approach for obtaining an approximate solution 
of the thermal boundary layer was described in reference 1. 
It is still simpler than the approach by Kroujiline, since in 
this method it is not necessary to calculate the flow boundary 
layer prior to the determination of the thermal boundary 
layer. This method uses exact solutions of the boundary- 
layer e([uations which are known for a special type of pressure 
variation in the flow which is encountered on wedge-shaped 
bodies. These velocity profiles and the temperature i)rofiles 
are used to approximate the actual profdes for arbitrary 
pressure variations. A differential equation is set up 
with the condition that the growth of the boundary layer 


at any place on the C 3 dinder with arbitraiy cross section 
be the same as for the wedge-type flow when the bound ar^^- 
layer thiclmess and the pressure gradient have the same 
values in both cases. When this idea is applied to the 
momentum thiclviiess of the boundarv la}"er, the resulting 
equation is identical with Karnian’s integrated momentum 
equation. When it is applied to the convection thickness of 
the thermal boundary layer, the resulting differential 
equation fulfills the requirement that the heat transferred 
from the surface to the fluid must be found again within the 
boundary layer (ref. 11). This method was compared with 
exact solutions and experimental values in reference 1 and in 
investigations performed at the University of California 
(ref. 12), and the agreement obtained was quite satisfactory. 
This method can easil}’- be extended to include effects of 
variations of the surface temperature and of internal heating 
(ref. 11); however, the corresponding exact solutions for the 
wedge-type flow must be known. Such solutions, which take 
into account the effects of a surface temperature variation 
and of internal heating, are presented in references 11 and 13. 
The method ma}" also be extended to the transpiration 
cooling of porous surfaces as soon as the corresponding exact 
solutions for this case are Imowii. A few of these solutions 
are presented in reference 14. However, too few solutions 
are given for use as a basis for the approximate method. 
Brown (ref. 15) has recentU made an extensive calculation 
to obtain exact solutions in transpiration-cooled porous 
surfaces of the wedge-flow type which include the effect- of 
pressure gradients and of variable property values. The 
results of this calculation are now being used to extend the 
method mentioned in the preceding paragraph to transpira- 
tion cooling. 

This paper deals with the method described as the first 
type of approach. The method has the advantage of being 
applicable to cases for which the corresponding wedge-type 
flow and heat transfer are not known. 

STATEMENT OF PROBLEM AND SIMPLIFYING 
ASSUMPTIONS 

This rei)ort is a contribution to the problem of determining 
the development of the thermal boimdary layer and the 
heat-transfer coeflicients on a body of arbitrary cross section 
with porous walls in a two-dimensional flow. Figure 1 
shows the cross section of a bo(U of this type. At Keynolds 
numbers that are sufliciently high, the flow around the body 
may be subdivided into the boundaiy-la^^er region, which 
surrounds the body with a very small thiclmess, and the 



Figure 1. — Sketch of turbine blade indicating notations used. 
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potential flow, which determines the pressure distribution 
around the body. The highest pressure on the body is found 
at the stagnation point. The pressure then decreases in 
flow direction along both sides of the body and usually 
increases again later. The pressure variation along the 
body determines the development of the flow boundary la^^er 
and also whether the boundary layer is laminar or turbulent. 
Usually, the laminar part is confined to a region near the nose 
of the body. The laminar boundaiy-layer region is investi- 
gated herein. The flow of coolant through the porous 
surface may be described by the velocity with which the 
coolant leaves the surface in a normal direction. The ques- 
tions which will be answered in this report are (1) what are 
the local heat-transfer coefficient and the surface tempera- 
ture on any point along the body for any prescribed dis- 
tribution of tlie coolant velocity, and (2) what distribution 
of the coolant velocity gives a desired distribution of the 
heat-transfer coefficient and of the surface temperature 
around the body. Usually, for example, a constant wall 
temperature is most desirable and the problem is to determine 
that distribution of the coolant velocity which results in a 
constant wall temperature. 

A number of simplifications must be introduced in order to 
keep the time required for the solution of a special problem 
within a tolerable amount. The following assumptions are 
made: The flow is two-dimensional and in steady state, 
internal frictional heating within the boundary layer can be 
neglected, and property values (density, viscosity, and heat 
conductivity) may be considered constant. The influence 
of temperature gradients along the surface of the body is 
neglected. Although this influence may be considerable 
(ref. 13), there is no quantitative information available for 
transpiration-cooled surfaces that would permit this effect 
to be taken into account. 

In applying the method to the determination of the tem- 
peratures of transpiration-cooled turl)ine blades, neglecting 
internal frictional heating should be admissible, since the 
temperature differences within the l)oundarv layer generated 
by the cooling process are considerably larger than those 
generated by aerodynamic heating. This will l)e shown in 
more detail later. On the othei* haiid, tlie large temperature 
variation in the boundary layer is connected with a con- 
siderable variation of the property values. The influence of 
this variation on the heat transfer may be approximately 
corrected by use of the results in reference 15. 

BOUNDARY-LAYER EQUATIONS 

In a coordinate system, showm in figure 1, where the :r-axis 
runs along the surface and the y-direction is normal to the 
surface and under the assumptions mentioned in the previous 
section, the differential equations describing the velocity and 
the temperature within the boundary layer arc 


(All symbols are defined in ajipendix A.) The equations 
have to fulfill the following boundary conditions: At the sur- 
face {y=0)y the velocity component u* parallel to the surface 
is 0, the velocity component v* normal to the surface lias a 
finite value Vu,*, and the temperature of the surface is inter- 
connected with the velocity by the over-all heat-transfer 
process, so that only one of the two values may be prescribed. 
At the outer edge of the boundary layer {y= oo), the u* com- 
])onent of tlie velocity transforms asymptotically into the 
stream velocity U* and the temperature transforms into the 
sti’eam temperature. Since only temperature difl'erences 
appear in equation (3), the temperature level does not enter 
into the problem. Therefore, all temperatures will be meas- 
ured from the teni])erature in the stream as reference tem- 
perature, and t will be interpreted as the temperature differ- 
ence fi*om this reference temperature. Consequently, outside 
of the ])oundary layer, ^ = 0. 

In a gas stream, it is advantageous to interpret t as the 
total temperature difference. In this case, equation (3) 
already includes the effect of the frictional heating for a gas 
with a Prandtl number of 1. Since for all gases, the Prandtl 
number does not deviate much from the value 1, equation 
(3) also gives a good appi-oximation to the real conditions for 
gases as long as the temperature differences impressed upon 
the boundary layer by a cooling jirocess are larger than the 
temperature differences created by internal friction. 

In order to reduce the numlier of parameters, the differen- 
tial equations will be made dimensionless. For this purpose, 
all lengths measured parallel to x are divided by i, the dis- 
tance between the stagnation point and the trailing edge of 
the body measured along the surface, and all velocities in 
this direction are divided by an upstream velocity Uq. All 
lengths and all velocities parallel to y are, in addition, multi- 
plied by the square root of Reynolds numbers Re^ based on 
the body length L and the upstream velocity C^o* 

2/=xV^ (4) 


The pressure gradient bp^jbx'^is impressed upon the boundary 
layer by the potential flow outside the boundary layer and 
can be replaced by the stream velocity gradient by use of 
Bernoulli’s equation: 


dp* 

dx* 




( 5 ) 


In this way, oquations (1) to (3) transform into 

bu , du TrbU , 
bx by bx by^ 


bu 

bx 


+ 



( 6 ) 

( 7 ) 


bt bt^ I bH 
^ bx ^ by Pr bi/ 


( 8 ) 


I 


* 


du* 

dy* 


dp* d^u* 
dx* ^ dy*^ 


du* dv* 
dx* dy* 


dt j dH 


( 1 ) 

( 2 ) 

( 3 ) 


The boundary conditions for these equations liecome 

7/=0 U = 0 V=Vu>{x) t = t^{x) (9) 

7/-^ GO u-^U ^-^0 (10) 

The equations are now integrated over y from y=0 to y= °o . 
The result of this integration is Karman’s integrated momen- 
tum equation 
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Tz ''»+S “ '> 

and the heat-flow equation 

L Jo (|^)„ 


d 

(lx , 


( 12 ) 


In order to simplify these ecpiations, the following eharae- 
teristic houndary-layer thieknesses are introduced: 

(1) The dis])laeement thickness of tlie flow l)oundary 
layer 




(13) 


(2) Tlie momentum oi* im])ulse thickness of the flow 
l)oundary layer 

(3) The convection thickness of the thei-inal boundary 
layer 


'•'Jo U 


(ly 


(I 5) 


where d is the ratio t/tu, with the limiting values 0=1 for 
y = 0 and 0=0 for Introducing these boundary- 

layer thicknesses and writing dUldx as transforms equa- 
tions (11) and (12) into 


0 


d 

dx 


(16) 

( 17 ) 


VELOCITY AND TEMPERATURE PROFILES 


In oi’der to obtain an a])])]‘oxiniate solution of the last two 
(Hjuations, a])])roximate ex])ressions for the velocity and 
temj)ei*ature ])rofdes will be inti’oduced. The accnirac.y of 
the ]-esults of this calculation will de])end on how well the 
actual pi’ofiles are a])])roximated by the assumed sha])es. 
The tem])erature ])rofiles as well as the velocity profiles are 
chosen as a one-i)aranieter family. Tlie i)arameter for each 
family is deterjnined in such a way that the assumed ])rofiles 
fulfdl the exact l)oundary-layer eciuations (6) to (8) at the 
wall surface: 




/^\ _i_ /£>'( 


(18) 

(19) 


In addition, (lu>, following boundary conditions will be 
fidfilled: 

?/=() M =0 0=1 ( 20 ) 


y- 


u^iJ 0^0 


( 21 ) 


Velocity profile. — For A>0, 0<^<1, 


u 

U 


=sin||+A(l-e"^^-sin||) 


Originally, polynomial expressions were used to a])proximate 
both the velocity and the temperature profdes. Schlichting, 
however, iiointed out that better approximations may be 
obtained liy ex])j*essions coni])osed of trigonometric functions 
(ref. 6). The following jirofiles are used in this report: 


and for y > 1 , 
0 


u - 3 - 

^=l-Ae ■ 


( 22 ) 


By use of this jirofile, the following expression is found from 
equation (18) for the shape parameter A: 


A=- 




9+ 


(:>-i). 


(23) 


Tbc vi'locity jirolile c((uation (22) lias ab-eady been used by 
Schlichting (ref. 6). For a solid flat plate (L"' = 0 and Vu=^), 
it approximates the Blasins velocity profile very well, and for 
V = 0 and t>„= — 3/5, transforms into the exact asymp- 
totic suction profile as calculated by Schlichting (ref. 6). 

For A<0, 0<f<l, 

0 

y \ K ' V { ^ y\ 


and for 7 > 1, 
0 


^=1 

u 


(24) 


The shape parameter in (bis case, as obtained from equation 
(IS), is 


vJ 


.■X- 5 

TT+R Vu-5 


(25) 


For the solid Hat ])late, this ])rotile gives the same expression 
as equation (22). It will be shown later that the profile 
a])i)roximatcs the sejiaration profile as calcidated by Ilartree 
(ref. 16) better than the usual polynomial expression. 
Sejiaration ocenrs at A= — 1 and the corresponding profile 
has the shajie 

(26) 


U • 2 ”■ ?/ 
2 6 


Temperature profile. — For 0< — <1, 

Ot 

t ^ . . TT y . . TT y / . . 7T y\ 


and for ^ > 1 , 

Ot 


0=0 


(27) 


The shai)c factor for this i)rofilc is obtained from equation 
(19) as follows: 
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1 


I 


^ _ —Pr vjt 
^ T-\-Pr Vu,dt 


(28) 


In this case, a negative shape factor is always connected 
with a positive value of Vw, wliereas for the velocity profile, 
the sign of the shape factor depends, in addition, on the 
pressure gradient dUldx. 

With these profiles, the different boundary-layer thick- 
nesses and the velocity gradient at the wall may be calculated. 




Y=0.1 366+0.03791 A-0.00786A2 




>-A>0 


5" 


|=0.1366-0.01456A-0.02618A^ ^A<0 


{ du\ T^U f . 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 


The nondimensional temperature gradient at the wall is 


TRANSFORMATION OF BOUNDARY-LAYER EQUATIONS 

FLOW BOUNDARY LAYER 

^Multiplication of equation (16) by dijU and a partial differ- 
entiation of the first term give 

The expression in the parentheses of tlie second term on the 
left-hand side and the term on the right-hand side are func- 
tions of the shape parameter A. Therefore, the momentum 
equation for the boundary layer may be written in its final 
form. 

|^=/2(A)-/.(A)C/'6i^+«,„5, (37) 


with the following expressions for the two functions and 
/ 2 , whicli are obtained from equations (29) to (34): 


,//A)=2+ 


0.3634-0.03005A 


0.1366+0.03791 A— 0.00786A2 
/2(A)=|^^+(3-5)Aj(0.1366+0.03791A-0.00786A2 
0.3634— 0.1366A 


M>0 


/i(A)=2- 


0.1366-0.01456A-0.02618A^ 


/2(A)=2 (1 +A) (0.1366-0.01456A-0.02618A^ 


VA<0 


(38) 

(39) 

(40) 

(41) 



-fi(l+A0 


(35) 


The corresponding expressions for the convection thickness 
of the thermal boundary layer are pj'esented in ap])ondix B. 


These functions are presented in figure 2. All curves have 
a break at A = 0 because different expressions approximate 
the velocity profiles for positive and negative A values. 
Equation (37) is a linear first-order differential equation 
from whicli the momentum thickness of the boundary layer 


.I9r 


. 10 - 


.17 - 


. 16 - 


h. 

s 




.( 5 - 


.14 - 


.13 - 


. 12 1 - 



Figurb 2. — Chart for determination of /i. A, and 5*75 used in flow boundary-layer differential equation. (An 18- by 10}^-in. working chart of this 

figure may be obtained upon request from NACA Headquarters.) 






6 


REPORT 1220 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



-1.0 -9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 


I .2 .3 .4 .5 


Figure 3. — Cliart for det(‘rmination of shaj^e parameter A for flow boundary-layer calculations. (A 22- by 15V2-in. working chart of this figure 

may be obtained upon request from NACA Headquarters.) 


can l)(‘ ohlahied ])V integral ion, as soon as tiie gradient L' 
of the stream velocity U and the i)oroiis How characterized 
1)V Vw fiiT known as functions of x. In order to make the 
calculations more convenient, the shape parameter A can be 
ex])ressed as a function of the two quantities U'bi and 


Vy,bi 


AX),A= 


J T/ ^ O 5 . bl 

bi Vjcbi ~ 


(42) 


A<0,A 


v,A ^ 


h 

5 



[ ^ !■ 
d“ 2 


b 


(43) 


The ratio bi/b occurring in these equations is a function of A 
(see eqs. (30) and (33) and fig. 2). The functional relation 
for A is plotted in figure 3. By use of figures 2 and 3, the 
integration procedure for the differential equation (37) be- 
comes very simple. The step-by-step procedure for such a 
calculation is presented in appendix C. 

In order to start the calculation at the stagnation point, 
the boundary-layer tliickness bf at this location must l)c 
known. At the stagnation point, the stream velocity U is 
zero. Since, on a blunt nose, the increase of the boundary- 


layer thickness is never infinite, the term on the right-hand 
side of equation (37) has to l)e zero. This gives the equation 

foi^) — ^^A(A)f = 0 (44) 


However, figure 3 also applies to the conditions at the stag- 
nation point. From both relations, U^bi^j ^,/n and /2 were 



I' KJURE 4. — Chart for determination of U'5i^ at stagnation point for 
flow boundary-layer calculations. 
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each obtained as a function of v^bi hy a trial-and-error 
process. These values are presented in table I. The con- 
nection between U'bi^ A, and Vw^i is also shown in figure 3 
as the stagnation line. Table I cannot be used immediately 
to start the calculation because the value v^bi is not known. 

However, the value which contains only known 

U bi U 

values, is also a function of U'bi^ and is plotted in figure 4. 
P>om this figure, U'bi^ can be determined for a given vJ!^' 
and bi can then be obtained from this value. 


THERMAL BOUNDARY LAYER 

With the help of equations (28) and (35), the heat-flow 
equation (17) is transformed into 


_1 X . , 7t 1 / PrKvy.b \ Ubt,c dt^ 

Vwdx^ 2 PrKvu,d\ ir+PrKvwb/ Vwtw d^ 


(45) 


where K denotes the ratio btjb. The last term on the right- 
hand side accounts for a variation of the temperature along 
the surface. A difficulty arises in connection witli this term. 
It is knowm from the results of references 11 and 13 that a 
temperature gradient dtjdx also has a i)ronounced effect on 
the shape of the temperature profile. Since no such effect 
was included in the assumed profile (eq. (27)), the signifi- 
cance of the last term is doubtful. It was therefore neg- 
lected herein, restricting this report to cases where the 
variation of the surface temperature is kept small either by 
internal conduction within the wall or by proper choice of 
Vw The investigation of the influence of large surface tem- 
perature gradients will be left to future work. With this 
simplification, equation (45) can be \n*itten in its final form 

Tx =^-MPrKvJ) (46) 


where /3 is a function of the product Pr as presented in 
figure 5. For a solid wall (Vu:=0), the equation simplifies to 


d 

dx 


KPrb 


(47) 


The shape parameter A^ may be written in the form 

_ —Pr Kvrcb 
^ 7T+Pr Kvw^ 


(48) 


This equation, together with the expressions for bt,d^ 
(appendix B), determines a functional relation between 
/i, bt,d^^ A, and Pr which is presented m figiu*e 6. The 
figure presents the ratios bt,d^ for each of three values of the 
shape parameter A (—1, 0, and +1). For intermediate 
values of A, linear interpolation in the range A= — 1 to 
A = 0, or A=0 to A=l, at a constant value of K may be 
used with good accuracy. The heat-flow equation (46) 
is again a linear first-order diff'erential equation from which 
the thermal convection thickness is obtained by integration 
when the stream velocity (7, the flow through the porous 
surface and the Prandtl number Pr are prescribed, and 
the boundary layer thickness b and the shape factor of the 
flow boundary layer A are knowm from a preceding solution 


of equation (37). Tlie step-by-step procedure for such a 
calculation is explained in appendix C. 

In order to start the calculation, the convection thickness 
bt^c niiist be known at the stagnation point. A partial 
differentiation of the left-hand side of equation (46) gives 

At the stagnation point, U is zero, and the boundary-layer 
increase dbt,ddx is not infinite on a blunt-nosed body. This 
results in the equation 

r ID 

or 

Kv^S) (49) 

VwO 0 



Pr 


Figure 5. — Chart for determination of fz used in thermal boundary- 
layer differential equation. 


428231—57 2 
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Figure 6. — Chart for determination of K for thermal boundary-layer calculations. (A 17- by 20-in. working chart of this figure may be obtained 

upon request from NACA Headquarters.) 
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This is a relation between the two unknown ratios bt,d^ and 
K. A second relation is given by figure 6. From both, the 
values and K can be determined by a trial-and-error 
procedure. The values are presented in table II and 
figure 7. 

When the thermal convection boundary-layer thickness 
is known, the local heat-transfer coefficient follows by a 
simple calculation. The equation which defines the local 
heat-transfer coefficient h is 




(50) 


This equation gives, for the local Nusselt num])er based on 
the length Z, 




(51) 


The introduction of equation (17) for the temperature 
gradient leads to tlie expression 


Nu _ 


=Pr [i dx Pr (f,- 


1) (52) 


Without porous cooling, equation (52) simplifies to 
Nu 7T 1 


(53) 


In many cases, all the heat transferred from the outside 
flow to the wall surface is picked up by the cooling fluid on 
its flow through the porous wall. In this case, the wall 
surface temperature may be calculated by the equation 


hi IP — pCpViQ (tp td) 


(54) 


where Z is the temperature with whicli the coolant enters 
the porous wall. The ratio of the difference in gas temper- 
ature minus wall surface temperature to the difference in 
gas temperature minus coolant temperature is given by the 
expression 




1 


h 


1 

>3 


(55) 


pCpVy; 


COMPARISON OF ASSUMED VELOCITY AND TEMPERATURE 
PROFILES WITH EXACT SOLUTIONS 

The accuracy of the method depends on how well the 
assumed profiles approximate the actual ones. It is there- 
fore necessary to check the accuracy of this approximation 
with the results of exact calculations to the extent that these 
are available. Such a comparison will be made in this section. 

VELOCITY PROFILES 

For the solid surface (z?^=0), the comparison can be made 
with exact solutions, which were obtained by Hartree (ref. 
16), for wedge- type flow for which the stream velocity is a 
power function of the distance from the stagnation point 
{U=UiX'^). A set of velocity profiles taken from reference 
16 is presented as a family of dashed lines in figures, where 
the ratio of the velocity u in the boundary layer to the stream 
velocity U outside of the boundary layer is plotted over the 
dimensionless distance 


2 /* y lU 

V2— iS V px* V2— isV X 





Figure 7. — Chart for determination of and fz at stagnation point for thermal boundary-layer calculations. (A 17- by 9-in. working chart 

of this figure may be obtained upon request from NACA Headquarters.) 
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Figure 8 . — Comparison of approximate velocity profiles for solid 
surface with Ilartree solutions (ref. 16). Vw, 0. 


and tlie shape parameter is defined by the equation 

^ 171 / 

i3 =— 7 • The expression for the dimensionless displace- 

m + 1 

d lU 

ment thickness is 2 i=- — — -y/ — Therefore, 

2 — jS V ^ 




2-/3 


U 

X 


U (3 U 

A dilferentiation of U=UiX”^ gives U'=rn —=- -- — 

® X 2-/3 X 

Therefore, 

(56) 


The value Zi as a function of ^ is presented in reference 1. 
Since is a function of A, according to equations (42) 

and (43), equation (56) presents a relation between /3 and A 
from which A may be calculated for any value of Tlie 
shape of the approximate velocity profile is determined by 
A (eqs. (22) and (24)). In this way, the solid profiles of 
figure (S were determined. They are siqierimposed on the 
exact profiles in such a way tliat the displacement thickness 
is the same in both cases. Figure 8 shows tbat generally 
the agreement between the exact profiles and the ajiproxi- 
mation used in this report is satisfactory. Only for the 
separation profile (/3= — 0.1988) are the deviations larger, 
but even for this separation profile the approximation by 
equation (24) is better than the usual four-term polynomial 
apjiroximation, which is shown as a dash-dot line. 

For a ])orous wall, some exact solutions are contained in 
reference 14. The notations in this reference are the same 
as the ones mentioned in the i)revious i)aragraph in connec- 
tion with the Ilartree solutions. In that report the porous 
flow velocity Vu? is cliaracteriziHl by 

From this and the previous expression for Zt 

^2—p X Zi=vJ>i (57) 

The value of Zi is presented in reference 14 as a function of 
X for two values of (3 (0 and 1). Therefore, ecpiation (57) 
presents a relation between X and A for a certain value of p. 
Equation (57) therefore connects the parameter X, which 
determines the shape of the exact velocity profile, with the 



z 


Figure 9. — Compari.son of approximate velocity i)rofiles near stag- 
nation point with Schlichting solutions (ref. 14). /3, 1. 



Figure 10. — Comparison of approximate velocity profiles for fiat 
plate with Schlichting solutions (ref. 14). 0. 


form parameter A determining the shape of the apiiroxima- 
tion. Figures 9 and 10 show a comparison between the 
exact profiles and the approximations for the neighborhood 
of the stagnation point (/3=1) and for a fiat plate (0=0). 
It may be seen that again this agreement is quite satisfactory 
for /3=1. For 0=0 the a])proximation is not as good for 
profiles with a distinct S shape. For the largest coolant 
flow (characterized by the highest value of X or the smallest 
value of A) the deviation is consideralile, altliough the agree- 
ment is somewhat better than that with the four-term poly- 
nomial a])proximation indicated by a dash-dot line. There- 
fore, the method presented in this report sliould not be used 
for excessive coolant-flow rates. 

TEMPERATURE PROFILE 

For a solid wall (r^j^O), the approximate expression 
(eq. (27)) for the temperature profile has a unique form, 
independent of the pressure gradient along the surface. 
For a gas with a Prandtl number of 1 and constant pressure 
along the siudace, it is known that the temperature profile 
is similar to the Blasius velocity jirofile. This similarity 
also holds for the approximation in this report. In addi- 
tion, it is shown in reference 1 that the shape of the actual 
temperature profile is influenced only to a minor degree by 
a pressure gradient in the flow. Therefore, these approxi- 
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Figure 11. — Comparison of approximate temperature profiles for 
flow at stagnation point with exact solutions (ref. 14), 




(a) Blade I. Ratio of chord length to distance between consecutive 

blades c/s, 1.408. 

(b) Blade II. Ratio of chord length to distance between consecutive 

blades c/s, 1.136. 

Figure 12. — Turbine-blade profiles used for calculated examples. 

Illations should be quite good for a solid surface. For a 
porous wall, with a constant pressure along its surface, 
and for a fluid with a Prandtl number of 1, the velocity and 
temperature profiles are again similar. Therefoi*e, figui-e 10 
also shows the degree of approximation for the temperature 
profile. 

For the flow in the neighborhood of a stagnation point 
and a fluid with a Prandtl number of 1, some temperature 
profiles were calculated in reference 14. Figure 11 presents 
the exact temperature profiles and the approximations in a 
way that shows the thickness of the temperature boundary 
layer, which is defined in the same way as the dis])lacenient 
thickness, to be the same for the accurate solution and the 
corresponding approximation. Figure 11 shows tliat the 
agreement is satisfactory as long as the coolant flow char- 
acterized by the value X or is not too large. The afiproxi- 




(a) Blade 1. 

(b) Blade ll. 


Figure 13. — Velocity distribution around surface of blade. 

Illation for the highest coolant flow rate (X = 3.191) is not 
shown since it is obvious that the exact curve cannot be 
approximated by the trigonometric functions from which 
the ap])roxiniation is composed. Large coolant flows in 
the ranges where the agreement ceases to exist will probably 
not be used because the gain in cooling effect for a given 
increase in coolant flow is too small in this range. 

NUMERICAL CALCULATIONS FOR TURBINE BLADES 

The method developed in the jirevioiis sections was used 
to calculate the local distributions of the heat-transfer 
coefficients for two turbine blade shapes for which the velocity 
distributions were known. Tlie two blades are shown in 
figure 12. The line at the nose of the blade indicates the 
circumferential direction and the arrow, the direction of the 
upstream velocity. Figure 13 shows the velocity distri- 
butions around the blade circumferences. The velocity 
distributions around the blades were calculated using the 
method described in reference 17. 

The method outlined in this report can be used to calcu- 
late the development of the laminar-flow boundary layer 
in the downstream direction up to the point of separation. 
In reality, at the Kcynolds numbers encountered in gas- 
turbine engines, the boundary layer usually becomes turbu- 
lent before it reaches this point. The transition point to 
turbulence is determined by the pressure gradients along 
the blade surface, by the boundary-layer thickness, by the 
temperature distribution within the boundary layer, and, in 
addition, by the curvature of the blade surface and by the 
stream turbulence (ref. 18). Little quantitative knowledge 
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(a) Blade I. 

(b) l^lade II. 

Figure 14 . — Growth of flow impulse and thermal convection boundary- 
layer thicknesses for blade. Pr, 0.7. 

exists on the iiiniieiiee of all these parameters. However, 
there aj*e indications that on turbine blades the ])oint of trans- 
ition is near the point at which the maximum value of the 
velocity is reached. A calculation which deals with laminar 
])oundary layers is therefore useful for the region near the 
nose of a turbine lilade. The points of maximum velocity 
are indicated in figures 14 to 16 by small dots. 

The first part of the calculation deals with the prolilem 
in which a uniform coolant flow rate Vw is ])rescribed along 
the l)lade periphery. The develo])inent of the flow and of 
the tlu'rmal houndaiT layers was calculated in the way 
described in a])peiidix C. The residts are ])resented in 
figure 14. A])parently l)oth boundary layers start out with 
a finite A^alue at the stagnation point and incj*ease in thick- 
ness in the downstream direction on l)oth the suction and the 
l)ressure siu’faces. The boundary-layer thicknesses at any 
])lace along the blade increase with increasing coolant flow 
rate. 

In figure 16 the local heat-transfer coeflicients along the 
blade suj*faccs tliat are olitained from the thermal boundary- 
layei- thickness are shown. Tlie figures show that very 
high local values are encountered at the stagnation point 
and that the values decrease consideralily with increasing 
distance fj*om this point. The ajiplication of transpii*ation 
cooling deci'cases the heat -transfer coefficients effectively 
from the values obtained on a solid blade surface (rt 4 ,= 0). 
The decrease is not as pronounced at the stagnation point 
itself as on the sides of the lilade. The value of the velocity 
ratio r«;*/Fo* may be obtained from the iiarameter 
presented in the diagrams hy dividing it by the square root 
of the upstream Keynolds number. For turbine blades 
this Reynolds number is near the value 10^, so that a velocity 

V * 

ratio 0.003 1 6 corresponds to a value ofru?=l. Values 


for this velocity ratio which are considered in practical 
api)lications are in the neighborhood of 0.01. A comparison 
of figures 15 (a) and (b) shows that generally the values of 

Xul-^jReo at the stagnation point decrease considerably 
with an increasing ratio of the radius of curvature at the 
blade nose to the blade chord. The values of N'uf^Reo are 
determined by the velocity gradient dUldx^ which also 
depends mainly on the dimensionless radius of curvature. 
Increasing the radius of curvature at the blade nose is 
therefore an effective means of decreasing the heat-transfer 
coefficient in this region. Downstream of the point of 
boundary-layer transition to turbulence, the heat-transfer 
coefficients will increase. To give some indication of the 
magnitude of the values that may be expected in the turbu- 
lent region, heat-transfer coefficients were calculated at the 
location a; = 0.5 under the assumption that these values are 
the same as on a flat ])late at the same distance x from the 
leading edge, with a stream velocity equal to the local value 
found at the blade at o; = 0.5, and with an upstream Reo 
equal to 10^ Formulas derived by Rannie and Friedman 
(ref. 18) were used for these calculations. The short horizon- 
tal lines in figure 15 indicate these values; the dashed part 
of the curve Vw=0 in figure 15 (a) shows the probable heat- 
transfer values in the transition and turbulent regions. 

The suj-face temperature of the l)lade is determined by 
the heat-transfer coefficients. Equation (55) gives the 
relation between l)oth values when heat conduction in 



(a) Blade I. 

(b) Blade II. 

Figure 15. — Variation of local heat-transfer coefficient near stagnation 
point for blade. Pr. 0.7 
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(a) Blade I. 

(b) Blade II. 

Figure 16. — Surface temperature distribution for blade. 0.7. 

the blade wall and lieat radiation can be neglected. Figure 16 
shows the values 1— $ determined from equation (55). 
The quantity 1 — is the ratio of the difference in blade 
surface temperature minus coolant temperature to the 
difference in effective gas temperature minus coolant tem- 
perature. The higher heat-transfer coefficients near the 
leading edge create higli blade temperatures at that location. 
Heat conduction within the blade wall tends to reduce 
tliese high local temperatures. This reduction is assisted by 
the fact that the blade temiieratures are especially low 
in the laminar region on both sides of the leading edge. 
Therefore, it may be that the blade-nose temperature is 
reduced to a value not higher than the blade-wall temperature 
in the turbulent part, especially for blades with a blunt nose. 
On the other liand, it has to be kept in mind that the heat- 
conductivity values for porous materials are lower tlian for 
solid walls (ref. 18). It has already been mentioned that 
temperature gradients along the blade surface influence the 
local heat-transfer coefficients somewhat; however, not 
enough information is available to account for this effect. 


A surface temperature that decreases in How du-ection 
tends to increase the heat-transfer coefficients so that the 
decrease in surface temperature on both sides of the stagna- 
tion point may be slightly less than those shown in figure 16 
where heat conduction within the blade walls is excluded. 

A second set of calculations was made for blade I to deter- 
mine that distribution of the local coolant flow rate v^o along 
the blade surface which results in a constant blade tempera- 
ture. The procedure for such a calculation is also described 
in appendix C. The local coolant flow rates which are 
necessary in order to keep the temperature I’atio ^ (eq. (55)) 
at the value 0.7 are shown in figure 17. The flow boimdaiy 
layer and the thermal boundary layer determined in this 
way are iiresented in figure 18. A coin])arison of this figure 
with figure 14(a) indicates that the development of the 
boundary layers is considerably different for both cases. 
The boundary-layer growth is smaller for the case of con- 
stant wall temperature. Figure 17 shows that the highest 
local flow rates are necessary near the stagnation point in 
order to keep the wall temperature down at that place. The 



Figure 17. — Variation of coolant flow required to maintain constant 
blade temperature for blade I. <l>, 0.7; Pr, 0.7. 



layer thicknesses for constant blade temperature for blade I. 
0.7; Pr, 0.7. 
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ma^nitiulo of the coolant flow rate at the stag;nation point 
is proportional to the square root of the velocity gradient 
(lUldx, which is itself determined mainly by the value of the 
radius of curvature at the blade nose. Tlie smaller this 
radius, the larger the velocity gradient and the local coolant- 
flow rat(\ The flow rates decrease on both the jiressure and 
the suction sindaces in the downstream direction. Down- 
stream of the transition point the boundary layer will be 
turbulent. It is known that the cooling eirectiveness of a 
turbuhuit boundary layer is less than that of a laminar 
boundary layer. Correspondingly higher values of the cool- 
ant flow rate are necessary in the turbulent region in o]*der 
to keej) the blade tenqierature constant. In order to give 
some indication of the magnitude of the coolant flow rate 
necessary for the turlndent region, the values Vu: which result 
in tlie same value (^ = 0.7) as in the laminar region were 
calculated at the location j==0.5. For this calculation the 
formulas jiresented in reference 18 were used and it was again 
assumed that the local <I> value on the blade surface is the 
same as on a flat plate at the same distance from the leading 
(Hlge and with a velocity equal to the local stream velocity 
V and an upstream Reynolds number 7?Co equal to 10k The 
flow i*ates obtained in this way are indicated by short 
horizontal lines in figure 17. The values in the tuiRulent 
region along the blade can be expected not to deviate much 
from tlu'se values. Apparently a considerably higher coolant 
flow rate is necessaiy in the tuiRulent portion of the blade 
surface than in most of the laminar iiortion. In order to 
simplify manufacturing problems, turbine blades may be 
produced in such a way that the required coolant flow rate 
has a constant value along the ])ressure surface and a differ- 
ent constant value along the suction surface. This then re- 
sults in lower local blade temperatures in the region of laminar 
boundaiy layers; the heat conduction into this cooler part 
of the blade may reduce the local teniiieratures at the stag- 
nation ])oint to values equal to or lower than the temjiera- 
tures in the turbulent region, especially when the radius of 
curvature at the blade nose is not very small. For blades 
with a very i)ointed nose it may be necessary to have a 
coolant flow rate at the stagnation point appi*oximately 
twice as large as along the sides in order to reduce the tein- 
])erature at that point. 


CONCLUSIONS 

A method which permits the approximate calculation of 
local heat-transfer coefficients and surface temperatures in 
the laminar flow region around cylinders of arbitrary cross 
section with transpiration-coolod walls was developed. Veloc- 
ity and teni])erature profiles in the boundary layer were 
a])])roximated by trigonometric expressions. The method 
was applied to determine local heat-transfer coefficients and 
surface temperatures in the laminar region of two transpi- 
ration-cooled turbine blades for a given coolant flow rate. 
Coolant-flow distributions necessary for maintaining uni- 
form blade temjieratures were also determined. The follow- 
ing conclusions are made: 

1. For small coolant flowrates, the assumed trigonometric 
a])])roximations for the velocity and temperature profiles 
within the boundary layer agree well with exact solutions. 
For high coolant flow rates (which are outside the range of 
])ractical interest), the agreement liecomes poor. It is, how- 
ever, better than the four-term i)olynomial approximation 
formerly used. 

2. Transpiration cooling results in a considerable reduction 
in heat-transfer coefficients for turbine blades in the laminar 
region. 

3. The surface temperatures of a blade with negligible 
heat conduction are highest at the leading edge, lowest in 
the rest of the laminar region, and have intermediate values 
in the turbulent region. 

4. The following variation of the coolant flow rate along 
tlie blade surface is necessary to keep the blade temperature 
constant. Highest local coolant flow rates are required at 
the stagnation point in order to keep the blade surface tem- 
perature down at that point. The coolant flow rates de- 
crease very markedly in the downstream direction on both 
the suction and pressure surfaces in the laminar region. 
However, they increase again when the turbulent region is 
reached. 

5. In order to maintain the same surface temperature on 
both sides of the blade, different coolant flow rates should be 
used for the suction and pressure siu’faces. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 22, 1951 


APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 
c chord length of turbine blade, ft 

Cp specific heat at constant pressure, Btu/(lb)(°F) 
fi function of shape factor of velocity profile (see eqs. 

(38) and (40)) 

/2 function of shape factor of velocity profile (see eqs. 

(39) and (41)) 

/g function of Pr, K, and h (see eq. (46)) 
h heat-transfer coefficient, Btu/(sec) (sq ft) (°F) 

K ratio of thermal to flow boundary-layer thickness, 
dtid 

k thermal conductivity, Btu/(sec) (ft) (°F) 

L length of either side of surface of profile measured 
from stagnation point, ft 

m exponent (termed Euler number) in U=UiX^ 

Nu Nusselt number, hL/k 

Pr Fraud tl number, Cpiijk 

p* pressure along body, Ib/sq ft 
Reo Reynolds number, UqLIv 

s distance between blades, ft 

t temperature difference between local and frec-stream 

values, °F 

U nondimensional stream velocity, U*IUq 

U* stream velocity, ft/sec 

Uq upstream velocity, ft/sec 

Ui constant in wedge-flow velocity, U=Uix”^ 
u nondimensional velocity component along surface, 
uVUo 

u* component of velocity along surface, ft/sec 
V nondimensional velocity component normal to sur- 

/ — 

face; yj-^Rcq 
Uq 

z;* velocity component normal to surface, ft/sec 

X nondimensional distance along surface, x*jL 

a;* distance from stagnation point along surface of 
profile, ft 

y nondimensional distance normal to surface, y-yP^o 

y* distance normal to surface, ft 




fi 

A 


5 


nondimensional distance used in wedge-type flow, 

y lU 

■yj2 — fi\ X 

nondimensional displacement tliickness used in 
wedgo-type 

shape parameter for wedge-type flow, 2m/(m+l) 
shape parameter of velocity profiles (see eqs. (23) 
and (25)) 

shape parameter of temperature profile (see eq. (28) 
or (48)) 

flow boundary-layer thickness 

flow boundarv-Iayer displacement thickness, 

/;(-&> 

flow boundajy-layer momentum or im]:>ulso thick- 


ness 


’Jo t/0 t/)*' 


5, 

5/. e 


thermal boundary-layer thickness 

thermal ])oundary-layer convection thickness^ 




nondimensional temperature ratio, tjtw 

shape parameter for exact solutions for porous walls, 


-r-r 


e 

X 

' X 

\'V 

fjL viscosity, ll)/(ft) (sec) 

V kinematic viscosity, m/P; sq ft/sec 

p density, Ib/cu ft 

4> nondimensional temperature ratio, ttjtc 

wSu])scripts: 

c coolant (at blade surface) 

w wall 

0 upstream condition 

wSuperscripts: 

* dimensional 

' difl'erentiated with respect to x 


APPENDIX B 

DETERMINATION OF THERMAL BOUNDARY-LAYER CONVECTION THICKNESS 

The determination of the thermal boundary-layer convection thickness 8t,o which is similar to the impulse thickness of 
the flow boundary layer, results from the evaluation of equation (15) 




u 


dy 


(15) 


when the assumed temperature and velocity profiles are inserted in the integrand. A single temperature profile was assumed; 

y 

it is given by equation (27) and is for 0 < ^ < 1 


o . . 7T 7/ , • 7T y • 7T 7/\ 


and for > 1 

Ot 


( 27 ) 


0=0 
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Distinct velocity profiles were assumed, depending on the sign of the shape parameter A. Moreover, these profiles were fiir- 

y y 

filer subdivided into different expressions depending on^ < 1 or ^^1- These assumed velocity profiles are given by equations 

(22) and (24); they are for A>0, 0<7<1. 

0 




and for^ > 1 




u 


(22) 


for A<0, 0<|<1 


u 

Tl 


. 7T 7/ , ^ . 7T 7/ / ^ . 7T 7/\ 

=sin-|+Asin-^(^l-sin-^) 


and for % >l 


^=1 

U 


(24) 


It is at once obvious that, for the evaluation of bt,c, four distinct cases must be considered. 

The upper limit of the integral in the evaluation of 5^,^ luay be restricted to the value bt by virtue of the assumed tempera- 
ture profile. If bt<bj the value of bt,c can be determined from a single int(‘gral with 0 and bt as limits, because the integrand 
vanishes for y>bt. On the other hand, if two integrals are required for the evaluation of bt,c because of the assumed 

velocity profiles. These integrals have limits 0 and 5, and b and bt. 

Results of the evaluation of bt,c for the four cases follow: 

Case I— A<0 and bt<b: 


5 5 A 5 !;)][*“ i 'i '^^ i I (' i !)] 

Integration and simplification lead to the following result: 

“y (1 +A) (l— cos /iC)— ^ (a"— 1 sin 7 tA')+(1 +A+A,+AA,) 


dy 


(A+AA,)[ 


K- 

ir(2/iC-l) + 7r(2K+l) T(2K-\f^^ 




.COS'] 


+ 


^‘^[t(2-] 


1 


(A,+AA,)| 


2) w(2K+]) 

“Y)“ir^+7v) (* +2“)]+ 


AA, 


[A sin ;tA-8^(/^^ ,) sin x(A> 1 Sin 


87r(/iC+l)‘ 

Case II— A<0 and 5,:>6: 

J„-JJ[l-sm:2-A,si,.:S(l-sm|2)][sin|f+Asin||(l-sini|)]*+ 

X*'[‘ i s]-^' i s] (' 1 0] <» ‘‘’■I 


Integration and simplification in this case lead to 
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^-(1 +A) 1 (1 +A,) cos 


/K \.K . 


(1 +a+A,+AA,){;^j^ sin g sin g ^^]| + 


(A+AA()| 


_ p_ 

T^{2K^\yir{2K-\) 


r7r(2i^-l)1 fxi 

\_ 2K J 7t(2/i:+1)'='’®L 


fA -LAA ^ 

, 1 

1 r 

>(2-A)-| 

1 fTT 

(A‘+AA,)^^(2+A) 

' 5t(2-A) 

7r(2-/7)‘^°®L 

. 2/7 J 

7r(2+A)^®®L 


77 

r-(/7^-i_ 

77 

^.r+A+i)-i) 

A 4 

85r(77— 1) 

L K J 

87t(A+1)"‘"L K J( 


^2K+l)-|_2g . ) 

2K \ -K 2K]^ 

Ik 


Case III — A>0 and 5^ <5; 

i i i 0] l f -r^-sin 5 1)] c/2/ 

Integration and simplification for tliis case lead to 

^c=(l_A) l-il-A) I cos "f +A/ir+§ c-3--§+(-l +A-A.+AA,) 


. ir{K-\) 

.1)^1“ K 


K 


sin J (77+1)] 


+ 


r 2 A+ 2 K 12/73 r 1 , 1 

' ^ -J+(a,-AA,) 


(A+AAO |_(3ei73+;r3) T (36 ^H,t3) 

1 /"i ,77\ 2 ^/7"|, 

"2 J+^^‘ b 


1 


K 


6(9/73++)^6(9/73++) 


K) w{2—K) 


-ZK I 3773 

® ■'■(9773+^2) e 


'] 


Case IV — A>0 and 5, >5: 

i i I I 0] I f+^ (l-e“^-sin 1 1)] C/ 2 /+ 

I i I i 0] 

Integration and simplification for this case result in 


, 2 , 
d 


. 2A , 1 /I I ^ \ 2 /l 7 T , A /K \ K . 7 t\ 

+Y+A — 1-(1+A,) — cos 2/^^+3 c +A, t~2 — 2+^ K ^ 


(-1+A-A,+AA,) 


r 2 K^ 


r 77 . 1 

(77-1)1 

1 • 1 

(A+1)1 

b(77-l)®“'| 

L 2 ‘ 'A J 

1 +77+1) 

L 2 A J 


(A+AA,) 


12 A3 


2/7 2/7 TV 

L(36A3++) ^ + ^‘^®®2/7 (36A3++) 


'] 


1 


1 


AA, 


ri A . TT 
[2~^ A 


(2+A) ^ir(2-A) ir{2-K) 
x3 , 3A3 


COS 


7t(2 — A) 

~2K 


+ 

1 


CK COS 


^(2+A)-l 


-3K_ 


6(9A3++) 


^(2+77) 2 A J 

] 


+ 


A 6(9A3++)^(9/73++) 

APPENDIX C 

PROCEDURE FOR CALCULATION OF FLOW AND THERMAL BOUNDARY LAYERS 

In tliis appendix tlie procedure will be explained by wliicli the flow and thermal boundary layers, the heat-transfer 
coefficients, and the temperature distribution around a body of given shape may be obtained. As outlined previously, the 
boundary conditions on the surface of the body may be prescribed in two different ways, either by prescribing the coolant 
velocity or by prescribing the temperature around the siu'face. The calculation procedure becomes simpler in the 
first case and this will be considered first. In addition to the value v^*, the distribution of the stream velocity U* just outside 
the boundary layer around the body must bo known. This velocity distribution may be obtained by any of the known calcula- 
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tion procedures for frictionless flow on bodies without flow separation or it has to be determined experimentally when flow 
separation occurs. From the velocity U'^{x) the gradient can be found. 

The values mentioned have to be transformed into the dimensionless quantities U{x)^ U'{x), and by use of equations (4). 

FLOW BOUNDARY LAYER 

The aim of this calculation is to determine the momentum thickness 8i as a function of the distance x from tlie stagnation 
point as measured along the surface of the body. 

As a first step, the value 8i has to be found at the stagnation point. For this purpose, determine vJIU' , and read the 
value U'8i^ from figure 4, thereby determining the momentum thickness 8i. The intersection of with tlie stagnation line 
in figure 3 determines A, the shape parameter for the velocity profile. The velocity ])rofile itself may be determined from 
equation (22) or equation (24) if it is desired. The required parameters may also be taken from table I. 

The momentum thickness along the profile surface is found from equation (37). This first-order linear differential 
equation may be solved by any of the known procedures. The method of isoclines was used for the solution of the numerical 
examples described previously, and its use will be outlined here. Equation (37) determines the direction of the tangents to 
the different 8i^ curves which satisfy this equation. The task is to find the one curve which contains the 8i value calculated 
previously for the stagnation point. Figure 19 shows the 8i^ values as ordinate and the distance x from the stagnation point 
as abs(ussa. The dii-ections of the tangents may l)e obtained from equation (37) for any point in this figure characterized by 
a pair of values x and 8i^ and may be inserted on the ordinates through a chosen sequence of distances Xi, X 2 , x^ . . . along 
the abscissa for several 8i^ values. The calculation proceeds for x and a chosen 8i in the following way: determine U'8i^ and 
Vj^8i and read from figure 3 the value A and from figure 2,/i and/ 2 . Calculate from equation (37), d(8i^)ldx. Now the direc- 
tion of the tangent may be inserted in figure 19. The same calculation is re])eatcd for other 8i^ values and the corresponding 
directions of the tangents are inserted in figure 19. Tlie same calculation is repeated for the distances X 2 , x^ . . . If the 
distance between these values is chosen small enough, a curve tliat starts out with the predetermined value 8f^ at the stagna 
tion point and is tangent to the straight lines inserted in figure 19 can be drami in the figure with good accuracy. The calcu- 
lation may be shortened considerably when the curve is inserted step by step after the tangents have been calculated for any 
value X, Then the correct value 8f^ for the following distance x may l)e guessed from the shape of the curve uj) to that point 
and the tangents need to be calculated only in the neighborhood of this value. 

The values which will be needed for the calculation may be arrangcnl in a table such as the following: 


X 

u 

U' 

vlp 

1 

1 

A 

8 

0 

.Tl 

X 2 








The momentum thickness 8i in this table results from the jireceding calculation. The shape parameter A was found in the 
course of this calcidation and the value 8 may be determined from figure 2. 

THERMAL BOUNDARY LAYER 

The calculation now has to l)e restricted to a fluid witli a certain Prandtl number. The value which characterizes the 
thermal boundary layer and will be determined in this section is the convection boundary-layer thickness 8t,c- To find this- 



X\ X 2 ^3 ^4 

Figure 19. — Isocline solution of flow boundary-layer equation. 
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value at the stagnation point, determine from the preceding table, read from figure 7 or table II btjb and /a. This deter- 
mines bt,c- Equation (52) determines the local Nusselt number Nu, and for the case in wliich the heat-balance equation (54) 
is satisfied, equation (55) determines also the temperature ratio 4>. The shape factor determining the shape of the tempera- 
ture profile according to equation (27) may be determined from figure 5 and equation (48). 

The convection boundary-layer thickness along the surface of the profile is determined by equation (46). This equation 
may be solved in the same way as the corresponding flow equation. Figure 20 indicates the procedure. The product Ubt,c 
is plotted over the distance x from the stagnation point. The short straight lines in the figure again indicate the direction 
of the tangents. It is known from exact solutions that the gradient of any boundary-layer thickness at the stagnation point 
is zero and that the boundary-layer thickness itself has a finite value. Thus, the curve Ubt,c in figure 20 starts out with the 
value zero and with the inclination U^bt,c for x=0. In order to find the direction of the tangents for any value a;, for example 
Xi, assume a value Ubt,c, calculate bt,c and btjb, read the value K from figure 6, determine the product Pr Vu>bK, read /s from 

figure 5, and calculate from the heat-flow equation (46) the gradient ^ (Ubt,c)- Insert the corresponding tangent into figure 20 

and repeat the calculation for other values Ubt,c- Now the curve which begins with the value zero for a;=0 and with 

the corresponding inclination, may be extended to Xi, The calculation is now repeated for the next distance X 2 , the curve is 
extended, and this procedure is repeated until the curve is known for the whole length x. The result of this calculation gives 
bt,c and /a. These values may be inserted in a table such as the following: 


X 

&C.C 

/a 

Nu 

4> 

Xi 

X 2 

3^3 






The ratio NulVP^oj which characterizes the local heat-transfer coefiScient, is found from equation (52). For the special case 
that the heat-balance equation (54) is fulfilled, equation (55) determines the temperature ratio <l>. 



Figure 20. — Isocline solution of thermal boundary-layer equation. 
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The calculation is more tedious wlieu the surface temperature instead of the flow velocity is prescribed along the blade. 
When the lieat-balance equation (54) is satisfied, the temperature ratio may be determined as a function of x. Now the j 
inomeiitum equation (37) and the heat-flow equation (46) have to be solved simultaneously by a trial-and-error procedure. I 
A value Vv> is assumed. The calculation procedure described jireviously is carried out to determine a temperature ratio 4> I 
corresponding to the assumed The calculation has to be repeated until the determined 4> value matches the prescribed 
value. This calculation has to be carried out stepwise for the distances x^^ 0 ^ 3 , .. . ! 

In some cases heat may be transferred to the surface of the profile by radiation and may flow along the wall of the body i 
by conduction. Then an equation which takes these processes into account has to re])lace the equation (54). The determi- | 
nation of the coolant flow velocity tv which results in a prescribed blade surface temperature, may be determined in the same j 
way as was just described. The procedure, however, becomes rather tedious in this case. I 
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TABLE I— STAGNATION POINT VALUES FOR FLOW BOUNDARY LAYER 


Vn}8i 

A 

U'8i^ 

ii 

8 

8 

/. 

h 

Vyj8 

- 0. 4999 

1. 0000 

0. 0001 

0. 1667 

0. 3333 

4. 0000 

0. 5000 

-2. 9988 

4 

0. 78 

. 0008 

. 1614 

. 3399 

4. 1059 

. 4336 

-2. 4783 

3 

. 6065 

. 01935 

. 1567 

. 3452 

4. 2029 

. 3820 

-1. 9145 

2 

. 486 

. 0345 

. 1532 

. 3488 

4. 2768 

. 3471 

-1. 3055 

1 

. 411 

. 05225 

. 1509 

. 3510 

4. 3260 

. 3257 

-0. 6627 

0 

. 359 

. 0715 

. 1492 

. 3526 

4. 3633 

. 3110 

0 

. 05 

. 329 

. 0812 

. 1484 

. 3533 

4. 3800 

. 3050 

. 3380 

. 1 

. 3075 

. 09025 

. 1475 

. 3541 

4. 4007 

. 2966 

. 6780 

. 2 

. 279 

. 1105 

. 1466 

. 3550 

4. 4216 

. 2888 

1. 3643 

. 3 

. 242 

. 130 

. 1453 

. 3561 

4. 4508 

. 2785 

2. 0647 

. 4 

. 225 

. 151 

. 1447 

. 3566 

4. 4644 

. 2739 

2. 7643 

. 5 

. 208 

. 1718 

. 1441 

. 3571 

4. 4781 

. 2692 

2. 4698 

. 6 

. 1835 

. 19175 

. 1433 

. 3579 

4. 4976 

. 2627 

4. 1870 

. 7 

. 172 

. 213 

. 1429 

. 3582 

4. 5066 

. 2596 

4. 8985 

. 8 

. 161 

. 234 

. 1425 

. 3585 

4. 5158 

. 2566 

5. 6140 

. 9 

. 151 

. 255 

. 1421 

. 3588 

4. 5250 

. 2539 

6. 3336 

1. 0 

. 142 

. 276 

. 1418 

. 3591 

4. 5324 

. 2515 

7. 0522 

1. 1 

. 135 

. 2975 

. 1416 

. 3593 

4. 5374 

. 2497 

7. 7684 

1. 2 

. 128 

. 3195 

. 1413 

. 3595 

4. 5442 

. 2478 

8. 4926 

1. 3 

. 123 

. 341 

. 1411 

. 3597 

4. 5492 

. 2464 

9. 2133 

1. 4 

. 114 

. 362 

. 1408 

. 3599 

4. 5561 

. 2449 

9. 9432 

1. 5 

. 113 

. 384 

. 1408 

. 3600 

4. 5568 

. 2439 

10. 6534 

1. 6 

. 107 

. 405 

. 1406 

. 3602 

4. 5618 

. 2424 

11. 3798 

1. 7 

. 101 

. 426 

. 1403 

. 3604 

4. 5688 

. 2406 

12. 1169 


TABLE II— STAGNATION POINT VALUES FOR THERMAL BOUNDARY LAYER 



II 

p 

Pr=l 

II 

to 

II 

o 

K 

8 

h 

K 

^t.c 

8 

h 


8 

h 

K 

8 

h 

0 

1. 80 

0. 3888 

CO 

1. 56 

0. 3140 

CO 

1. 2 

1 

0. 2042 

00 

0. 66 

0. 0742 

i 

0. 05 

1. 76 

. 3949 

4. 40 

1. 54 

. 3271 

1 3. 610 

1. 215 

. 2286 

, 2. 5448 

. 78 

. 1207 

1.317 

. 1 

1. 72 

. 4142 

2. 555 

1. 55 

. 3477 

' 2. 113 

1. 26 

I . 2616 

1. 607 

. 94 

. 1773 

: 1. 081 

. 2 

1. 74 

. 4291 

1. 629 

1. 58 

. 3808 

1. 408 

1. 35 

. 3172 

1. 198 




. 3 

1. 74 

. 4518 

1. 336 

1. 61 

. 4124 

i 1. 233 

1. 43 

. 3664 

1. 091 




. 4 

1. 74 

. 4690 

1. 228 

1. 64 

. 4371 

1. 142 

1. 50 

. 4027 

1. 053 



! 

. 5 

1. 74 

. 4869 

1. 158 

1. 66 

. 4607 

1. 096 

1. 57 

. 4338 

1. 0325 




. 6 

1. 75 

. 4995 

1. 116 

I. 68 

. 4783 

1. 068 







. 7 

1. 78 

. 5101 

1. 087 

1. 70 

. 4933 ! 

1. 053 







. 8 

1. 79 

. 5208 

1. 068 

1. 72 

. 5066 

1. 040 





! 


. 9 

1. 80 

. 5295 

1. 056 

1. 74 

. 5174 

1. 032 







1. 0 

1. 80 

. 5371 

1. 048 

1. 76 

. 5266 

1. 025 







1. 1 

1. 84 

. 5444 

1. 038 










1. 2 

1. 84 

. 5488 

1. 037 








1 

1 


1. 3 

1. 84 

. 5524 

1. 028 










1. 4 

1. 85 

. 5569 

1. 024 










1. 5 

1. 85 

. 5608 

1. 021 








1 


1. 6 

1. 85 

. 5667 

1. 018 








1 


1. 7 

1 

1. 85 

. 5697 

1. 017 
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